
fer$Treatment[fer$Treatment=="Anoxic to oxic"&fer$Sample_day<=13]<-"Anoxic"
fer$Treatment[fer$Treatment=="Oxic to anoxic"&fer$Sample_day<=13]<-"Oxic"
f = fer%>%
filter(!is.na(Treatment))%>%
ggplot(aes(x = Sample_day, y = NPOC, color = Treatment))+
geom_point()+
ylab("NPOC (mg/L)")+
geom_vline(xintercept = 13.5)+
scale_color_manual(values = c("#175676","#D58936","#A44200","#4BA3C3")) +
facet_grid(Status~filt)+
theme_bw()
fer%>%
filter(!is.na(Treatment))%>%
ggplot(aes(x = Sample_day, y = TNb, color = Treatment))+
geom_point()+
ylab("Total bound nitrogen (mg/L)")+
xlab("Day of experiment")+
geom_vline(xintercept = 13.5)+
scale_color_manual(values = c("#175676","#D58936","#A44200","#4BA3C3")) +
facet_grid(Status~filt)+
theme_bw()
library(plotly)
ggplotly(f)
fer%>%
filter(Sample_day == 9,
filt == "Filtered",
Treatment == "Oxic")
fer
f

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkocmVhZHhsKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKCnNlcDA3IDwtIHJlYWQuY3N2KCIuLi9kYXRhL2NsZWFuX2RhdGEvbnBvY183c2VwMjEuY3N2IiklPiUKICBtdXRhdGUoRGF5ID0gYXMuRGF0ZSgiMjAyMS0wOS0wNyIpKQpzZXAwOCA8LSByZWFkLmNzdigiLi4vZGF0YS9jbGVhbl9kYXRhL25wb2NfOHNlcDIxLmNzdiIpJT4lCiAgbXV0YXRlKERheSA9IGFzLkRhdGUoIjIwMjEtMDktMDgiKSkKc2VwMDkgPC0gcmVhZC5jc3YoIi4uL2RhdGEvY2xlYW5fZGF0YS9ucG9jXzlzZXAyMS5jc3YiKSU+JQogIG11dGF0ZShEYXkgPSBhcy5EYXRlKCIyMDIxLTA5LTA5IikpCnNlcDE0IDwtIHJlYWQuY3N2KCIuLi9kYXRhL2NsZWFuX2RhdGEvbnBvY18xNHNlcDIxLmNzdiIpJT4lCiAgbXV0YXRlKERheSA9IGFzLkRhdGUoIjIwMjEtMDktMTQiKSkKc2VwMTUgPC0gcmVhZC5jc3YoIi4uL2RhdGEvY2xlYW5fZGF0YS9ucG9jXzE1c2VwMjEuY3N2IiklPiUKICBtdXRhdGUoRGF5ID0gYXMuRGF0ZSgiMjAyMS0wOS0xNSIpKQpzZXAxNiA8LSByZWFkLmNzdigiLi4vZGF0YS9jbGVhbl9kYXRhL25wb2NfMTZzZXAyMS5jc3YiKSU+JQogIG11dGF0ZShEYXkgPSBhcy5EYXRlKCIyMDIxLTA5LTE2IikpCnNlcDE3IDwtIHJlYWQuY3N2KCIuLi9kYXRhL2NsZWFuX2RhdGEvbnBvY18xN3NlcDIxLmNzdiIpJT4lCiAgbXV0YXRlKERheSA9IGFzLkRhdGUoIjIwMjEtMDktMTciKSkKc2VwMjAgPC0gcmVhZC5jc3YoIi4uL2RhdGEvY2xlYW5fZGF0YS9ucG9jXzIwc2VwMjEuY3N2IiklPiUKICBtdXRhdGUoRGF5ID0gYXMuRGF0ZSgiMjAyMS0wOS0yMCIpKQpzZXAyMSA8LSByZWFkLmNzdigiLi4vZGF0YS9jbGVhbl9kYXRhL25wb2NfMjFzZXAyMS5jc3YiKSU+JQogIG11dGF0ZShEYXkgPSBhcy5EYXRlKCIyMDIxLTA5LTIxIikpCnNlcDIzIDwtIHJlYWQuY3N2KCIuLi9kYXRhL2NsZWFuX2RhdGEvbnBvY18yM3NlcDIxLmNzdiIpJT4lCiAgbXV0YXRlKERheSA9IGFzLkRhdGUoIjIwMjEtMDktMjMiKSkKc2VwMjQgPC0gcmVhZC5jc3YoIi4uL2RhdGEvY2xlYW5fZGF0YS9ucG9jXzI0c2VwMjEuY3N2IiklPiUKICBtdXRhdGUoRGF5ID0gYXMuRGF0ZSgiMjAyMS0wOS0yNCIpKQoKZnVsbCA9IHNlcDA3JT4lCiAgZnVsbF9qb2luKHNlcDA4KSU+JQogIGZ1bGxfam9pbihzZXAwOSklPiUKICBmdWxsX2pvaW4oc2VwMTQpJT4lCiAgZnVsbF9qb2luKHNlcDE1KSU+JQogIGZ1bGxfam9pbihzZXAxNiklPiUKICBmdWxsX2pvaW4oc2VwMTcpJT4lCiAgZnVsbF9qb2luKHNlcDIwKSU+JQogIGZ1bGxfam9pbihzZXAyMSklPiUKICBmdWxsX2pvaW4oc2VwMjMpJT4lCiAgZnVsbF9qb2luKHNlcDI0KSU+JQogIGZpbHRlcighZ3JlcGwoInNwa3xkdXAiLFNhbXBsZSkpJT4lCiAgbXV0YXRlKFNhbXBsZSA9IHRvbG93ZXIoZ3N1YigiICIsIiIsU2FtcGxlKSkpCgpmdWxsJFNhbXBsZV9kYXkgPC0gYXMubnVtZXJpYyhzdWIoIlthLXpdK1swLTldK2QiLCIiLHN1YigiKFthLXpdKF8uKykqJCkiLCIiLGZ1bGwkU2FtcGxlKSkpCmZ1bGwkUmVwIDwtIGFzLm51bWVyaWMoc3Vic3RyKGZ1bGwkU2FtcGxlLDQsNCkpCmZ1bGwkZmlsdCA8LSBzdWIoIlthLXpdK1swLTldK1thLXpdK1swLTldKyIsIiIsc3ViKCJfLisiLCIiLGZ1bGwkU2FtcGxlKSkKZnVsbCRUcmVhdG1lbnQgPC0gc3Vic3RyKGZ1bGwkU2FtcGxlLDEsMikKZnVsbF9jb25jID0gZnVsbAoKZnVsbF9jb25jJFRyZWF0bWVudFtmdWxsX2NvbmMkVHJlYXRtZW50PT0iYWEiXTwtIkFub3hpYyIKZnVsbF9jb25jJFRyZWF0bWVudFtmdWxsX2NvbmMkVHJlYXRtZW50PT0ib28iXTwtIk94aWMiCmZ1bGxfY29uYyRUcmVhdG1lbnRbZnVsbF9jb25jJFRyZWF0bWVudD09Im9hIl08LSJPeGljIHRvIGFub3hpYyIKZnVsbF9jb25jJFRyZWF0bWVudFtmdWxsX2NvbmMkVHJlYXRtZW50PT0iYW8iXTwtIkFub3hpYyB0byBveGljIgoKZmVyID0gZnVsbF9jb25jCmZlciRmaWx0W2ZlciRmaWx0PT0icyJdPC0iRmlsdGVyZWQiCmZlciRmaWx0W2ZlciRmaWx0PT0idCJdPC0iVW5maWx0ZXJlZCIKCmZlciRTdGF0dXMgPC0gIk94aWMiCmZlciRTdGF0dXNbZmVyJFRyZWF0bWVudCAlaW4lYygiQW5veGljIHRvIG94aWMiLCJBbm94aWMiKV08LSAiQW5veGljIgoKd3JpdGUuY3N2KGZlciwiLi4vZGF0YS9wcm9jZXNzZWRfZGF0YS9UT0NfZm9ybWF0dGVkLmNzdiIsIHJvdy5uYW1lcyA9IEYpCgpmZXIkVHJlYXRtZW50W2ZlciRUcmVhdG1lbnQ9PSJBbm94aWMgdG8gb3hpYyImZmVyJFNhbXBsZV9kYXk8PTEzXTwtIkFub3hpYyIKZmVyJFRyZWF0bWVudFtmZXIkVHJlYXRtZW50PT0iT3hpYyB0byBhbm94aWMiJmZlciRTYW1wbGVfZGF5PD0xM108LSJPeGljIgoKZiA9IGZlciU+JQogIGZpbHRlcighaXMubmEoVHJlYXRtZW50KSklPiUKICBnZ3Bsb3QoYWVzKHggPSBTYW1wbGVfZGF5LCB5ID0gTlBPQywgY29sb3IgPSBUcmVhdG1lbnQpKSsKICBnZW9tX3BvaW50KCkrCiAgeWxhYigiTlBPQyAobWcvTCkiKSsKICBnZW9tX3ZsaW5lKHhpbnRlcmNlcHQgPSAxMy41KSsKICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gYygiIzE3NTY3NiIsIiNENTg5MzYiLCIjQTQ0MjAwIiwiIzRCQTNDMyIpKSArCiAgZmFjZXRfZ3JpZChTdGF0dXN+ZmlsdCkrCiAgdGhlbWVfYncoKQoKZmVyJT4lCiAgZmlsdGVyKCFpcy5uYShUcmVhdG1lbnQpKSU+JQogIGdncGxvdChhZXMoeCA9IFNhbXBsZV9kYXksIHkgPSBUTmIsIGNvbG9yID0gVHJlYXRtZW50KSkrCiAgZ2VvbV9wb2ludCgpKwogIHlsYWIoIlRvdGFsIGJvdW5kIG5pdHJvZ2VuIChtZy9MKSIpKwogIHhsYWIoIkRheSBvZiBleHBlcmltZW50IikrCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMTMuNSkrCiAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9IGMoIiMxNzU2NzYiLCIjRDU4OTM2IiwiI0E0NDIwMCIsIiM0QkEzQzMiKSkgKwogIGZhY2V0X2dyaWQoU3RhdHVzfmZpbHQpKwogIHRoZW1lX2J3KCkKCmxpYnJhcnkocGxvdGx5KQpnZ3Bsb3RseShmKQpmZXIlPiUKICBmaWx0ZXIoU2FtcGxlX2RheSA9PSA5LAogICAgICAgICBmaWx0ID09ICJGaWx0ZXJlZCIsCiAgICAgICAgIFRyZWF0bWVudCA9PSAiT3hpYyIpCgpmZXIKZgpgYGAKCg==